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O I The evolution of scalar perturbations is studied for 2-component (non-relativistic matter 

O I and dark energy) cosmological models at the linear and non-linear stages. The dark energy 
is assumed to be the scalar field with either classical or tachyonic Lagrangian and constant 
equation- of- state parameter w. The fields and potentials were reconstructed for the set of cos- 
p |. mological parameters derived from observations. The comparison of the calculated within these 
Q i models and observational large-scale structure characteristics is made. It is shown that for 
^ ' w = const such analysis can't remove the existing degeneracy of the dark energy models. 

^ ; Introduction 

> : 

0> ! The observations of the last decade surely confirm the acceleration of the cosmological ex- 
^ I pansion. The explanation of this fact needs the assumption that the main part - approxi- 
I mately 70% - of the energy density of the Universe belongs to the mysterious repulsive com- 
G^ ■ ponent called "dark energy". The simplest model describing satisfactory almost the whole 
^ ■ set of the experimental data is ACDM-one. Here dark energy is identified with the A-term 
(3 ! in the Einstein equations. However, in this case there are several interpretational problems, 
^ I which suggest that another solution should be found. The most popular alternative approaches 
I are quintessential scalar fields, i.e. scalar fields with the equation-of-state (EoS) parameter 
j_] ■ — 1 < Wde = Pde/ pde < —1/3. The simplcst physically-motivated Lagrangians are the classical 
. i and tachyonic ones. The first of them is the simple generalization of the non-relativistic par- 
ticle Lagrangian to the field while the second (called also the Dirac-Born-Infeld one) - of the 
relativistic particle one P, [HI ETJ ESI El]- The Lagrangian of classical field has the canonical 
kinetic term, the Lagrangian of tachyon field has the non-canonical one. 

As soon as the analysis of dynamics of expansion of the Universe [32] doesn't allow us to 
choose the most preferable by the observational data model of scalar field dark energy, here we 
focus on study of the evolution of scalar perturbations and the large-scale structure formation 
in the Universe filled only with the non-relativistic matter and either classical or tachyonic field 
minimally coupled to it. It should be noted that the behavior of perturbations has already been 
studied for different classical scalar fields more widely [5l[6l[34| than for tachyonic ones [Tl [TO]. 
The parametrizations of scalar fields, their impact on formation of the large-scale structure of 
the Universe as well as on cosmic microwave background anisotropy are widely discussed in the 
literature (see, for example, P, [131 [HI EH [30] ^ind citing therein). In this paper we analyse the 
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models with reconstructed for Wde = const potentials of the classical and tachyonic scalar fields 
and compare the obtained results to the ACDM-ones. 



1 Cosmological background 

We consider the homogeneous and isotropic flat Universe with metric of 4-space 

ds^ = Qijdx^dx^ = a?{r]) {drf — dapdx^dx^) , 

where the factor a(ri) is the scale factor, normalized to 1 at the current epoch rio, r] is conformal 
time {cdt = a{rj)drj). Here and below we put c = 1, so the time variable t = xq has the 
dimension of a length, and the latin indices i, j, ... run from to 3, the greek ones - over the 
spatial part of the metric: /i, ...=1, 2, 3. 

If the Universe is filled with non-relativistic matter (cold dark matter and baryons) and 
minimally coupled dark energy, the dynamics of its expansion is completely described by the 
Einstein equations 

7^,,-i^,,7^ = 87^G(7;f)+7:f)), (1) 

where TZij is the Ricci tensor and T^.j^\ ^ energy-momentum tensors of matter (m) and 
dark energy (de) . If these components interact only gravitationally then each of them satisfy 
differential energy-momentum conservation law separately: 

(here and below denotes the covariant derivative with respect to the coordinate x^). For 
the perfect fluid with density p[m,de) and pressure P{m,de), related by the equation of state 

P{m,de) = W(m,de)P{m,de), it givCS 

P{m,de) = -3-p(m,de)(l + W(^rn,de)) (3) 

a 

(here and below a dot over the variable denotes the derivative with respect to the conformal 
time: " '"= d/drj). The matter is considered to be non-relativistic, so w^n — and p^n — Pm oT^ 
(here and below "0" denotes the present values). 

We assume the dark energy to be a scalar held with either classical Lagrangian 

Caas^\(t>-^V-U{(t>) (4) 

or Dirac-Born-Infeld (tachyonic) one 



Ctach^-u{iUi-i4-^\ (5) 

where 0, ^ are the classical and tachyonic flelds respectively while W(0), U{^) are the held 
potentials deflning the models. We suppose also the background scalar flelds to be homogeneous, 
so their energy densities and pressures depend only on time: 

Pdas = + W(0) , Pclas = - W(0) , (6) 



Ptach = —F^ Ptach = I (7) 
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Then the conservation law gives the scalar field evolution equations 



) + 2ai^0 + a2^ = 0, (8) 
^-""^ ZaHi+^^^ = a, (9) 



1 



where H = d/a^ is the Hubble constant for any moment of the conformal time r]. 

We specify the model of each field using the EoS parameter Wde = Pde/ Pde- It is obvious 
that the scalar field evolution equations have the analytical solutions for w = const (here and 
below we omit index de denoting both - classical and tachyonic - scalar fields for Wde)- In this 
case another important thermodynamical parameter - the adiabatic sound speed = pde/ Pde 
- is equal to w. 

The analysis of the dynamics of the Universe expansion for the reconstructed fields with 
w = const was presented in |32]. It doesn't depend on the scalar field Lagrangian and - as a 
result - doesn't allow us to distinguish such models of scalar fields. So, in order to choose the 
most adequate to observations type of dark energy we should study at least the linear stage of 
the evolution of scalar perturbations. 



2 Evolution of scalar linear perturbations 

We derive the equations of evolution of scalar linear perturbations in dark energy - matter dom- 
inant era by varying of the Lagrange-Euler and Einstein equations in the conformal-Newtonian 
frame with space-time metric 

ds'^ = a^{ri)[{l + 2'<iJ (x, ri))dri^ - {l + 2^{x, ri))5a/3dx''dx% (10) 

where \&(x, 77) and ^(x,?]) are metric perturbations, which in the case of zero proper anisotropy 
of medium (as for dust matter and scalar fields) satisfy the condition \l/(x, i]) = — $(x, rj) exactly 
[UlTH]. In the theory of linear perturbations all spatially-dependent variables are usualy Fourier- 
transformed, so, all perturbations - of metric, fields, matter density and velocity - in equations 
are presented by their Fourier amplitudes: \1/(A;, r^)), (50(A;,?7), 5^{k,r]), 5^"^\k,r])^ V^'^'>{k,r]) 
etc., where k is the wave number. They are gauge-invariant - as it is particularly discussed in 
the original papers [4l[l5] and numerous reviews (see, for example, [8l[9l[25] and citing therein). 
The energy density and velocity perturbations of dark energy, 5^'^'^^ and V^^'^\ are connected 
with the perturbations of field variables 50, in the following way: 

.-.^(l + u.)(7-*+=ff). (11) 

yiclas) ^ 




w J u di 

yitach) ^ (14) 
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Other non- vanishing gauge-invariant perturbations of scalar field are isotropic pressure per- 
turbations 



and intrinsic entropy 



TT 



{das) _l+W I 6(t) 



w 



:i5) 



(16) 



Yi.de) ^ ^{de) _ ^^(de) 
^ W 



(17) 



The density perturbation of any component in the conformal-Newtonian gauge Ds = S, which 
is gauge-invariant variable, is related to the other gauge-invariant variables of density pertur- 



bations D and Dg as: 



D=D, + 3(l + .)(* + ^^)=B. + 3(l + .)^:^, 
where Dg, D, Dg and V correspond to either m- or de-component. 



(18) 



2.1 Evolution equations 

Evolution equations for scalar field perturbations S(f){k,ri) and 6^{k,ri) can be obtained either 
from Lagrange-Euler equation or from differential energy-momentum conservation law Tq.-^'^^^ — 
0: 



5(1) + 2aH5^ + 



d(p 



(19) 



2aH-,aH\l\ - --^ 



8i + 



Id^U fl dU 



U di^ \U di 



X 



5e - *e - 1 1 - 1 f ] ) + 2*^^ + = 0. (20) 



The linearised Einstein equations for gauge-invariant perturbations of metric and energy- 
momentum tensor components are 



^n^V^(™^+Pde(l + «^)V^('^)] = 0, 



+ ai/V"(™) - A;* = 0, 



+ = 0, 



d/'^'^ + 3(c^ - w)aEDf^ + k{\ + w)l/(''^) + 3aii"wr(''^) = 0, 



(21) 
(22) 
(23) 



V^"^^ + aE{\ - 3c2)\/('^^) - k{\ + 3c2)* - - J^r^'^^) = 0, (24) 



(25) 
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where 



w 



Yiclas) _ _ ^2) 



W 



Yitach) ^ _(y; + c2) 



A; 

k 



= (1 - c2)D(^''^^), (26) 
-{w + cl)D^"''''\ (27) 



In = const-case 1 — = 1 — w, w + c^ = 2w, hence the difference between equations for 
classical and tachyonic fields isn't big (for w close to —1 - as it follows from the observable 
data [T6j) and suggests the similarity of their solutions. 

So, in each case we have the system of 5 first-order ordinary differential equations for 5 
unknown functions \I/(A;, a), D^g^\k, a), V^"^\k, a), Df^\k, a) and V^'^'^\k, a) satisfying also the 
constraint equation: 



(28) 



2.2 Initial conditions 



Now we are going to specify the adiabatic initial conditions. The adiabaticity condition in 
two-component model gives Dg™-* = Df^^j (1 + w) [3, [8l [TS]. 

Since the density of the w = const-fields is negligible at the early epoch (a ^ 1), both 
our models are initially matter-dominated. It is known that in such case the growing mode 
corresponds to \E' = const. The field equations of motion for the reconstructed potentials are 
following: 



+ 



3w 



-3w 



9(1 -w) 




2 1 - + n^ea-^^ 



+ 



■3w] 



+ 



2 + w + w 



-3w 



1 ~ ^de + ^deO- 

rideil + w) 4a^' + 3(1 - u;)^ 



for the classical field and 



(29) 



6C - - + 3w; + 



3w 



-3ui 



2 1 - + ^dea-^"" 

^deO'~^^ 



6^' 



— w 



— 1 



5i 



w 



(1 - 3w)a^' - 



Hq\/\ — fide + ^^de«' 



-3ui 







(30) 



for the tachyonic one. Here and below a prime denotes the derivative with respect to the scale 
factor a and Qde = Pde/Pc, where pc = SHq/SttG. 
The condition = const for a <^ 1 gives: 




3 ifoVl — f^de 



de 



^ 3 



(31) 
(32) 
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Here T^'^^) = 0. 

Using these solutions and equations (fTTll-(fT4l). (fTSl) . (l28l) . one can find the initial values of 

Dl"'\k,a), V('")(A;,a), Df''\k,a), V^'''\k,a): 





2 k ^init 1 


(33) 


Jj(de) _ 

9 init 


-5(1 


(34) 


init 


2 ^ 


(35) 


£){m) _ 

9 init 




(36) 



which specify the growing mode of the adiabatic perturbations. 
2.3 Numerical analysis 

We have integrated numerically the systems of equations for dust matter and dark energy 
with w = const for the adiabatic initial conditions using the publicly available code DVERKlll. 
We used the set of cosmological parameters from http://lambda.gsfc.nasa.gov/product/map, 
assumed '^init = —1, o,init = 10~^° and integrated up to a = 1. The evolution of perturbations 
is scale dependent, so we performed calculations for k = 0.0001, 0.001, 0.01 and 0.1 Mpc~^. 
The models with the classical scalar field are denoted as QCDM, with tachyon - as TCDM. 
For comparison we also solve the evolution equations for the ACDM-model. 

As it can be seen in FiglH the simple conclusion, that the behavior of the scalar linear 
perturbations in the model with the tachyonic field with w = const should be similar to that 
in the model with the corresponding classical field, is valid. The curves for both fields almost 
overlap, so in this case it is not possible to choose the Lagrangian prefered by observations (see 
also [35]). 

In both models studied here the matter clusters while dark energy is smoothed out on sub- 
horizon scales. Generally, at present epoch the growth of the matter density perturbations is 
supressed and - unlike ACDM-case - such supression is scale dependent, however this depen- 
dence is very weak. The dark energy perturbations grow approximately up to the moment of 
the entering of particle horizon and start to decay after that (the density perturbation D^'^^^ - 
slowly) . 

Note that the perturbations in such scalar field models are insensitive to the initial condi- 
tions. Really, if we assume the dark energy to be initially homogeneous {6(f) = 6(j)' = and 
6^ = 6C,' = 0), the results of numerical integration will be the same as in adiabatic case (the 
similar conclusion was made in [H US]). 

The simplest test for identification of the source causing the accelerated expansion of the 
Universe can be based on the action of the studied fields on cosmic microwave background. 
Here the main attention should be paid to the temporal variation of the gravitational potential 
which causes the late-time integrated Sachs- Wolfe (ISW) effect. In Figl2]the evolution of \l/ is 
shown for both scalar field models and for ACDM-one for the scales of perturbations k = 0.0001, 
0.001, 0.01 and 0.1 Mpc~^. It can be easily seen that the scale dependence is weak and there 
is no substantial difference between classical and tachyonic dark energy. Such scalar fields are 
in many senses similar to the cosmological constant and their behavior is closer to that of the 
A-term for EoS parameter values closer to —1. The given dependence for ACDM-model doesn't 

^It is created by T.E. Hull, W.H.Enright, K.R. Jackson in 1976 and is available at 
http: / /www. cs.toronto.edu /NA / dverk.f.gz 
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Figure 1: Evolution of the density (top) and velocity (bottom) perturbations for the non- 
relativistic matter (left) and dark energy (right). In the left column the scales are k = 0.1, 
0.01, 0.001 and 0.0001 Mpc~^ from top to bottom. In the right column they are also k = 0.1, 
0.01, 0.001 and 0.0001 Mpc~^ from top to bottom for the density perturbations while for the 
velocity ones the curves correspond to A; = 0.001, 0.0001, 0.01 and 0.1 Mpc~^ from top to 
bottom at a ~ 1. The cosmological parameters are: fide = 0.722, w = —0.972, Qm = 0.278, 
h = 0.697. 



allow us to exclude this model using the observational data, because the difference between it 
and those in models with the scalar field dark energy is not substantial. 

It should be noted that neglecting of the dark energy perturbations leads to the "quasi- 
ACDM"-models, i.e. models in which the fields affect the growth of the matter perturbations 
only through the background. In these models the decay of the gravitational potential is scale 
independent and close to the small-scale one in the models with perturbed dark energy (in 
agreement with the results of [19j). 

Another possible test is based on the study of action of dark energy on the clustering proper- 
ties of dust matter. However, here we need the analysis of the evolution of scalar perturbations 
at the non-linear stage. 
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Figure 2: Evolution of the gravitational potential for the scales k = 0.0001, 0.001, 0.01 and 0.1 
Mpc~^ (from top to bottom). The curves for QCDM- (solid line) and TCDM- (dotted) models 
with the dark energy perturbations overlap. 

3 Spherical collapse in the models with homogeneous dark 
energy 

The simplest and most popular approach used in the study of the non-linear stage of the large- 
scale structure formation is the spherical collapse model. Within this framework we analyse 
the formation of the virialised haloes in the ACDM- and in the w = const QCDM- and TCDM- 
models with reconstructed potentials, discussed in the previous sections. 

The magnitudes of density perturbations of the classical and tachyonic scalar fields with scale 
less than the particle horizon are lower than corresponding magnitudes of the matter density 
ones by few orders and practically do not affect their growth. The amplitudes of matter density 
perturbations in the QCDM- and TCDM-models grow almost equally in cosmologies with the 
same parameters. They are also close to ones in the ACDM-models. So, in order to simplify 
the discussion of the non-linear evolution of scalar perturbations we assume the dark energy 
component to tend to homogeneity. Other reasons for homogeneous distribution of dark energy 
in the regions of matter inhomogeneties see, for example, in [22l [23] . Hence, the temporal 
dependence of the dark energy density is defined by the corresponding background equation. 
(Structure formation in inhomogeneous dark energy models has been analysed in [26].) 

The relative perturbation of mass of the dust component in the comoving volume v = 



47ri?V3 and metric ds^ = dt^ - M^{R)y\t, R)dR'^ - x^{t, R)R\cos'^ Odip^ + dO^) is following 




(37) 



where x(t, R) is the local scale factor derived from the Einstein equation Q 
the Einstein tensor and k is the Einstein constant) [T8l \T7\ : 




(here Qj is 




(38) 



(in this section 



a 



• 77 



d/Hodt). 
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For the A-term we should put Pde/ Pc = — while for the quintessential dark energy the 
relation is Pdel Pc = wflde- The local curvature parameter flf gives the amplitude of the initial 
perturbation: 6m(t) ~ |(^^a' — ^/)^m^o-(^) at a ^ 1. Since in the A-case Qf = ^f{R) for 
the dark energy we should put fi/ = Qf{t,R) [36l [17]. It means that here we need an addi- 
tional equation defining the evolution of the local curvature. However, for the homogeneous 
perturbations {-^^f = 0, x = x{t) = y{t)) using the combination of the Einstein equations 
Qq+Qi+ 2^1 = k(7^° + + 27^^) we obtain the motion equation without time-dependent 
curvature [36] : 

X 1 

- = (3pde + Pde + Pm) , (39) 

X Pc 

where the dust matter density is pm = pmX~^. Combining of the equation (l39l) and Friedmann 
equations for the homogeneous Universe allows us to find the evolution of the mass perturbation. 

In the analysis of halo formation the moment of reaching of the dynamical equilibrium is 
important. According to the virial theorem at this moment the kinetic energy becomes T^r = 
— ^Um,vir + UA^yir. lu the ACDM-modcl the energy conservation law T+Um + U\ = E is obtained 
by integration of the equation of motion (l38l) . multiplied by xx. From this follows: T = 
Um = —\^mX~^, Ua = — if^AX^ and E = ^Qf. In models with the A-term the total energy is 
constant in time and at the reaching of the dynamical equilibrium is equal E = \Um,vir+'^UA,vir- 
Alternatively, at the turnaround moment, when a; = T = 0, the total energy is E = Um,ta+U\^ta- 
From these equations we obtain finally: \Um,vir + '^UA,vir = Um,ta + f^A,ta- This identity, valid 
for the A-case, doesn't hold for the dark energy, since the temporal variation of the curvature 
in the perturbed region gives E{tmr) 7^ E(tta)- In other words, explicit temporal dependence of 
the dark energy density pde(t) leads to the explicit dependence of the potential energy of this 
component Udeit,x) on time and - hence - to the explicit temporal dependence of the total 
energy: E(t) = T{x)+Um{x) + Ude{t, x). It means that at different times we have different values 
of the total energy. For estimation of the moment of reaching of the dynamical equilibrium for 
the dark energy case we use the equations (3.13)-(3.17) from [17]. These equations describe 
the evolution of the spherically-symmetric perturbation with the arbitrary profile in the model 
with dark energy. Assuming there Qf = ftfit), x = x{t) = y{t) and V = V{t), we obtain the 
equations for the homogeneous spherical cloud. The additional condition of homogeneity of 
the dark energy (the equality of (l3|) and (3.17) from [17]) gives the expression for V. Using 
it together with (3.15) from [17] we obtain the equation describing the temporal variation of 
curvature: 



fi/ = 3 \{l + w)ndex\ (40) 

Combining of this equation with ( l38ll leads to the energy conservation equation Eif) = T{x) + 
Um{x) + Ude(t,x) in the following form: 

1 O 1 1 

2 2 x 2 2-^ ^ ' 

Really, it can be easily seen that differentiating (|4T1) with respect to time and using ((31), (l40l) we 
obtain (l39ll . Combining it with (|4T]1 we get (l38ll . Since E = ^fif, the equation (l40l) describes 
the temporal variation of the total energy. Using the virial theorem, for the moment of reaching 
of the dynamic equilibrium we write: E(tmr) = \Um,vir + '^Ude,vir- At the turnaround moment 
we have: E(tta) = Um,ta + Ude,ta- Taking into account that E(tmr) 7^ E{tta), we obtain finally: 

J^Um^vir '^U de,vir j-i/, \ (Um,ta ~l~ Ude,ta) • (^2) 

2 j^ytta) 
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This equality differs from used one in [20l[36] by factor E{tmr) / E{tta), which, however, is close 
to unity. 



0.008 




a(t)/a(U 



Figure 3: The temporal dependences of the local curvature parameter change 
(i7/(0) — Vtfit)) /Vtfifi) for 3 different f2/(0), which correspond to the redshifts of halo collapse 
Zcoi = 0, 1, 5. 

The mentioned above density contrast Ac is defined as the ratio of the density of viri- 
alised halo to the critical one at the expected moment of collapse: = Pvir/ Pc{tcoi) = 
QruXy^iHo/ H(tcoi)y ■ So, using the energy conservation law we find: 

'^^deitcol)xl„ + 2nf{t,ol)x^ir + = 0- (43) 

We choose the initial value of the curvature in the perturbed region ^f{0) by setting the 
moment of collapse tcoi- Unfortunately, in the scalar field plus CDM model the equation (l4Tll 
isn't symmetrical in time with respect to the moment of turnaround as it was for the ACDM- 
model. Hence the relation between the quantities tcoi and ^^/(O) must be established by the 
numerical integration of the equations (l40ll - (l4Tll . It is interesting to see how the curvature 
changes in the models with dark energy plus CDM, since in the ACDM-models it remains 
constant. In Figl3]we show such dependences for 3 different initial values of the local curvature 
parameter, which correspond to the redshifts of halo collapse Zcoi = 0, 1, 5. 

The calculations of Ac at the collapse moment Zcoi in the ACDM-model {fl\ = 0.721, 
nm = 0.279) and in scalar field plus CDM model (l^rfe = 0.722, w = -0.972, = 0.278) 
are presented in FigHl We see that the difference between the values of Ac in these models 
is insignificant for all Zcoi- It means that dynamics of cluster formation at all stages of their 
evolution - from linear stage through collapse up to virialization - is similar in the models with 
both reconstructed quintessential scalar fields if they are minimally coupled and have the same 
density (fide) and EoS (w) parameters. Practically, they are indistinguishable also from the 
best-fit ACDM-model, so, we conclude that these classes of dark energy models are degenerate 
with respect to their impact on dynamics of expansion of the Universe as well as formation of 
its large-scale structure. The other scalar field models (with special potentials, variable EoS 
parameter, non-minimal coupling, spatial inhomogeneity etc.) show similar but not so strongly 
degenerate impact on the linear and non-linear stages of structure formation (see, for example, 

[23 EH). 
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Figure 4: The density contrast Ac at the collapse moment Zcoi of spherical cloud in the 
Q(T)CDM- and ACDM-models. 

4 Theoretical predictions and observations 

In this section we compare predictions of the models with modern observational data on large- 
scale structure of the Universe. 

We have reconstructed the potentials of classical and tachyonic scalar fields for the model 
with parameters Vide = 0.722, w = -0.972, Vim = 0-278, = 0. 0467, h = 0.697, as = 0.799, 
Us = 0.962, taken hom\http:/ /lambda. gsfc.nasa.gov/product/ma]}i (see also [H]). For computa- 
tion of the CMB temperature and matter density fluctuations power spectra for TCDM-model 
we have modified the CMBFAST-4.5.1 code substituting the subroutines evaluating the classi- 
cal field perturbations by corresponding tachyonic ones. Here we used the homogeneous initial 
conditions for dark energy. 

The angular power spectra of CMB temperature fluctuations computed for QCDM- and 
TCDM-models are presented in the left panel of FiglHl There are overlaped curves for perturbed 
and unperturbed dark energy. For comparison the angular power spectrum of CMB temper- 
ature fluctuations for ACDM-model with parameters Vtde = 0.721, fi^ = 0.279, VLb = 0.0462, 
h = 0.701, (Ts = 0.817, Us = 0.96 [H] is computed and presented in the same figure. The 
corresponding five-year WMAP observational data [121 El] are shown there too. We have 
renormalized the computed power spectra by fitting them to all experimental points by x^- 
minimization procedure. The minimal in the ACDM-model equals 45.7, in the perturbed 
QCDM- and TCDM-models it equals 44.1 and in unperturbed ones 43.9. So, the difference 
between them for 43 experimental points is statistically insignificant. 

In the right panel of Figl5] the power spectra of matter density perturbations 
P{k) =< D^™) • > for the same models are shown. They are normalised to 5-year 

WMAP data in the way explained above. The power spetrum obtained from the analysis of 
the clustering of the luminous red galaxies in the Sloan Digital Sky Survey (SDSS LRG) |33| is 
shown by dots. The model spectra we fitted to observational one using scale independent bias 
parameter b [P°^^{k) = b'^P{k)) by the x^-minimization procedure. The differences between 
Xmin of 5 models are less than 1%, the best-fit bias parameter b equals 1.22 for ACDM-model 
and 1.24 for all scalar field plus CDM models analysed here. 

Thus, the ACDM-model as well as the perturbed and unperturbed QCDM- and TCDM- 
models can't be distinguished by current cosmological observational data. 
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Figure 5: The CMB temperature (left) and matter density (right) fluctuations power spectra for 
the models of Universe with dark energy: ACDM - solid line, QCDM - dashed, TCDM - dotted. 
The curves for perturbed and unperturbed dark energy (QCDM and TCDM) overlap with 
practically the same best-fit bias parameter (6 = 1.24). The observational CMB temperature 
and matter density fluctuations power spectra were obtained in the WMAP [24] and SDSS [33] 
projects. The amplitudes of the matter density fluctuations power spectra are normalised to 
WMAP 5-year data. 

Conclusion 

The evolution of the scalar perturbations is studied for the 2-component (dust matter and 
minimally coupled dark energy) cosmological models at the linear and non-linear stages. The 
dark energy component is assumed to be either classical or tachyonic scalar held with the 
potential reconstruted for the constant EoS parameter. 

The evolution of linear perturbations is similar for both types of Lagrangian. The small 
difference can be due to the generation of the intrinsic entropy of the fields, but in w = const- 
case it is almost the same for both types of dark energy as soon as the observational data prefer 
the values of the EoS parameter close to —1. 

The scalar fields studied here suppress the growth of matter density perturbations and the 
magnitude of gravitational potential. In these models - unlike ACDM ones - such suppression 
is weakly scale dependent and doesn't depend substantially on the Lagrangian. Such features 
can be used for calculations of the matter density power spectrum at different redshifts and 
of the power spectrum of CMB temperature fluctuations in the range of scale of the later 
integrated Sachs- Wolfe effect and - as a result - for interpretation of the observable data in 
order to identiflcate the nature of dark energy. However, the higher precision of the planned 
experiments could probably verify only whether the EoS parameter is —1 (ACDM-model) or not, 
but shouldn't remove the degeneracy of the dark energy models due to the type of Lagrangian 
for w = const. 
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